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Abstract 

Using an ab-initio molecular dynamics scheme (the Fireball scheme), we determined the equi- 
librium structure of intermediate size Si n Ge m (n + m = 71) nanoclusters with/without hydrogen 
passivation on the surface. Due to the strong surface distortion, defect states are found to permeate 
the energy gap of Si n Ge m clusters. However, the defect states are removed by adding H atoms on 
the surface of Si n Ge m clusters, and the gap opens up to a few eV, indicating a blueshift for photo- 
luminescence. It is also found that the radiative transition between the highest occupied molecular 
orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) states is enhanced by one 
to two orders of magnitude for Si n Ge m nanoclusters with respect to the corresponding pure Si 
clusters. This significant increase of the emission probability is attributed to the strong overlap of 
HOMO and LUMO wavefunctions that are centered mostly on the Ge atoms. 

PACS numbers: 73.22.-f, 61.46. +w, 71.15.Pd, 78.67.-n 
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I. INTRODUCTION 



The optical properties of bulk Si and Ge are rather mediocre because the light emission 
in the bulk Si and Ge is a phonon-assisted indirect process. Therefore, to improve on the 
light emission feature in Si-based materials is a challenge for both the technological and the 
fundamental research. 

The luminescence is a result of a significant overlap in electron and hole wave functions 
as the strength of the luminescence (i.e. the emission rate and quantum efficiency) depends 
on the extent of this overlap and the transition probability. Possible means to increase 
this overlap for the Si-based materials may be accomplished through, for example, alloying 
to change the band structure, impurities to produce the intermediate state through which 
the electron can recombine with the hole, or zone folding to yield the desired quasi-direct 
transition . However, the most important breakthrough in this issue is the observation of 
the visible photoluminescence (PL) from porous Si 0, 0, Q] and Si quantum dots P which 
opens new possibility for fabricating visible light emitting device from Si-based materials. 
Structural analysis of porous Si is quite difficult. But several measurements have confirmed 



that the principal feature of porous Si consists o 
enough to exhibit quantum confinement effects 
reported on Si nanowires 



^extremely fine structures which are small 
I. Various theoretical works have been 



11 



12, 



13. 



14. 



151 ] . They clarified 



|1CJ and Si clusters 

that the quantum confinement effects give rise to the change in the electronic structure and 
the optical properties and are the principle mechanism of the blueshift PL in porous Si and 
Si quantum dots. 

There were experimental reports indicating that Ge quantum dots embedded in Si02 
glassy matrices or in porous Si show a strong room t emp erature luminescence. 
Theoretical studies on the structure and stability of Ge clusters [18], the polarizabilities of 
small Ge clusters [n| , and the quantum confinement effect on excitons in Ge quantum dots 
|2fll | have also been reported. Since Ge has smaller electron and hole effective masses and a 
larger dielectric constant than the corresponding quantities for Si, the effective Bohr radius 



of the exciton in Ge is larger than that in Si, anc 



more pronounced in Ge than in Si 
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the quantum confinement effect appears 



2l|. These results suggest that Si n Ge r , 



nanoclusters could be possible candidates to be used as components for nanoscale functional 
optical devices. In order to understand the physics of Si n Ge m clusters, we performed an 
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ab initio molecular dynamics simulation for Si n Ge m clusters of an intermediate size, and 
systematically studied their electronic and optical properties. There is no doubt that the 
mismatch effect dominating the electronic and optical properties in Sii^Ge^ alloys may 
introduce interesting optical features in Si n Ge m clusters. But one has to keep in mind 
that the surface distortion associated with the stabilizing of the Si n Ge m clusters will also 
play an important role. Therefore the competition between the lattice mismatch and the 
surface distortion is the basic issue in our investigation. For this purpose, we studied Si n Ge m 
clusters of an intermediate size with/without the hydrogen passivation. By comparing our 
results between the two sets of Si n Ge m clusters, we found that the surface distortion plays an 
important role in the intermediate size of Si n Ge m clusters but the mismatch effect dominates 
when Si n Ge m clusters of the intermediate size are passivated by hydrogen atoms to eliminate 
the dangling bonds so as to lessen the effect of the surface distortion. We also found that 
the latter shows an enhancement of radiative transition and a blueshift in PL. 

It should be noted that our simulations of the Si n Ge m clusters without hydrogen pas- 
sivation will lead only to one of the more stable configurations among a large number of 
structural isomers. Hence its resulting structural and electronic properties may not ex- 
actly represent the corresponding properties of the true ground state configuration of the 
Si n Ge m clusters. However this caveat does not affect the main conclusions of our study, 
including the role played by the lattice mismatch and the surface relaxation in Si n Ge m clus- 
ters with/without hydrogen passivation, and the enhanced radiative transition in Si n Ge m 
clusters. 



II. METHOD 

We considered 71-atom Si n Ge m clusters, with m taking on the values 0, 18, 35, 53, and 
71, respectively. The initial configurations for the five clusters were generated randomly 
on a regular tetrahedral network. They were then relaxed by the ab initio molecular 



221. This 



dynamics scheme developed by Sankey and co-workers (the Fireball scheme) 
scheme is based on the density functional theory (DFT) in the local density approximation 
(LDA), where a local basis set is used to construct the Kohn-Sham orbitals. The basis 
functions are slightly excited pseudo-atomic orbitals (PAO). The Kohn-Sham orbitals are 
calculated self-consistently using the Hamann-Schliiter-Chiang pseudo-potential [23| and 
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the Ceper ley- Alder form of the exchange-correlation potential as parameterized by Perdew 



and Zunger 



24( | . We have tested this method in the study of the strain relaxation of 



Sii T Ge r alloys and the results are in good agreement with the experimental observations 



In our simulation, sp 3 -type PAOs were used with confinement radii of 5.0 ae, 5.2 cte, 
and 3.6 as for Si, Ge, and H atoms, respectively. A cubic cell with a lattice constant 50 
A was chosen as the unit cell. This cell size is sufficiently large to ensure that spurious 
interactions between clusters will vanish. To obtain the stable configuration of a Si n Ge m 
or a Si n Ge m H 84 cluster, molecular dynamics simulations were performed on the clusters, 
starting from the initial random configuration, at a temperature of 10 3 K for about 20 ps 
until the network was equilibrated. The network was then slowly cooled to 300 K for 5 ps, 
and the dynamical quenching was finally performed to fully relax the system to zero K. 
Charge transfer was calculated self-consistently in the simulations, which is important in 
modeling clusters. 



After performing the molecular dynamical simulations, we obtained stable configurations, 
five each for Si n Ge m and Si„Ge m H 8 4 clusters, with the ratio raj {n+m) = 0.0, 0.25, 0.51, 0.75, 
and 1.0, respectively. As an example, Fig. ^ (a) shows the stabilized configuration of the 
Si 36 Ge35 cluster. Apparently, there is strong surface distortion in this case and the stabilized 
structure is a compact network in a more oblate structure. This kind of structure has also 
been found to be more stable for the intermediate-size of Si clusters from other theoretical 
studies j^. It is found that the surface distortion leads to local bonding configurations 
with more than four bonds, in particular in the vicinity of the surface. Further structural 
analysis shows that the average bond lengths of Si-Si (bsisi), Si-Ge (bsiGe), and Ge-Ge (bceGe) 
of the Si n Ge m clusters are expanded relative to the bulk, as shown in Fig. |21 (a). But they 
still maintain the relationship of fesiSi < ^SiGe < &GeGe and are almost independent of the 
ratio of mj {n + m) . In addition, as shown in Fig. 121 (b), the average angles Q a ^i which 
shows the angle between Pa and fij atoms, are less than the tetrahedral angle, except ^siSiSi 
in SiigGe 5 3 cluster. The large fluctuation of the average angles 9 function of the 

ratio of m/(n + m) indicates that the clusters have lost their initial tetrahedral symmetry. 
This means that the surface distortion strongly dominates in the intermediate size Si n Ge m 
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clusters, whatever the ratio of Si/Ge. 

When the initial configurations were passivated with 84 H atoms to terminate the 
dangling bonds on the surface atoms, all the relaxed Si n Ge m H 84 clusters still maintain a 
tetrahedral symmetry in the interior, with only a minor distortion on the surface. Such 
features can be seen, for example, from the relaxed structure of Si36Ge3 5 H 84 shown in 
Fig. [T] (b). The average bond lengths bsisi, bsice, and bceGe of the cluster (see Fig. |21 (c)) 
and the average angles 9 a ^ (see Fig. |21 (d)) are quite close to the corresponding features 
in Sii^Ge^ alloys (comparing Figs. El (c) and (d) to Figs. 3 and 7 in ref. 2^|)- This is 
an indication that the surface distortion in hydrogen-passivated clusters is weak even for 
Si n Ge m H; clusters of an intermediate size. On the other hand, the mismatch effect becomes 
the dominating factor in hydrogenated Si n Ge m clusters, similar to the situation exhibited 
in Sii-^Gea; alloys j^ . 



III. ELECTRONIC STRUCTURE 

The LDA calculation is known to underestimate the energy gap in the single-particle 
energy spectrum of semiconductor clusters. Corrections to the underestimated energy gap 



can be calculated using the GW approach 



14 1271 28]. However, GW calculations are 



expensive, even for clusters of an intermediate size considered in this work (i.e., Si n Ge m Hg4 
clusters with n + m = 71). Our main interest in this study is to provide an understanding 
of the trend of change in the electronic structure and optical properties of Si n Ge m clusters 
of a given size as the concentration of the Ge component varies. The LDA approach is 
expected to be sufficiently accurate to provide a qualitative description of this trend of 
change. Therefore, we have calculated the electronic density of states (EDOS) of the relaxed 
Si n Ge m clusters, using the same DFT/LDA scheme as that in the structural determination. 

The EDOSs of Si n Ge m and Si n Ge m H 84 clusters at the five ratios of m/{n + m) described 
previously are shown in Fig. El- From the left panel of Fig. El many defect states associated 
with the surface distortion were seen in the HOMO-LUMO energy gap region of Si n Ge m 
clusters. These defect states were cleared up when the surface was passivated by H atoms 
(see the right panel of Fig. 3) because H atoms terminate the dangling bonds of the surface 
atoms, resulting in the clusters keeping the basic tetrahedral symmetry. Even though the 
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energy spectrum of Si n Ge m and that of Si n Ge m H 84 clusters are quite different, the general 
feature of the energy spectrum in both cases was found to be rather insensitive to the ratio 
of mj (n + m). 

Fig. E] (a) and (b) illustrate the HOMO-LUMO energy gap of Si n Ge m and Si n Ge m H 84 
clusters as a function of the ratio of m/(n + m), respectively. Because of the existence of 
the defect states in the energy gap region, the HOMO-LUMO gap of Si n Ge m clusters is 
small (only 0.4-0.5 eV) and does not show any dependence on the ratio of Si to Ge. But 
the HOMO-LUMO gap of the hydrogenated Si n Ge m clusters is opened up to several eV 
(3.5 - 4.0 eV) and shows a linear dependence on the ratio of Si to Ge which reflects the 



mismatch effect similar to that for the corresponding alloys [25]. Specifically, experimental 
observation [2^ and theoretical calculations jsoL l31 | indicate that the lattice mismatch 
associated with alloying in Sii_ x Ge x alloys induces a change of the conduction band, 
in particular, the lowest energy in the conduction band (LUMO) changes from the A 
point near X point for Si-rich alloys to the L point for Ge-rich alloys [3jJ, leading to a 
nearly linear decrease of the indirect band gap. In the case of Si n Ge m clusters passivated 
by hydrogen, the dominating factor is again the lattice mismatch. This similar effect 
is observed from the HOMO-LUMO gap as shown in Fig. |3Jd. The existence of such 
large energy gap for Si„Ge m clusters of this size (with n + m = 71) can be attributed 
to the effect of quantum confinement and is consistent with the scenario of a blueshift in PL. 



IV. OPTICAL PROPERTIES 

It is well known that an accurate calculation of the band gap is important for the determi- 
nation of the optical polarizability of semiconductors and insulators. Many efforts have been 
focused on the correction to the "too-small" band gap obtained by DFT/LDA calculations 
for the purpose of accurately predicting the optical response of these systems. Another 
key factor for determining the optical properties of these systems is a proper description 
of interacting electron- hole pairs (excitons). While calculations using the Green's function 
based on the GW approach have led to quite accurate prediction of the energy of low-lying 
excitations and the oscillator strength for small clusters [^, the application of this method 
to clusters of intermediate sizes is still computationally too expensive. Most recently, a 
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method based on the linear-response theory within the frame work of the time- dependent 
(TD) DFT/LDA had been applied to calculate the optical spectra of clusters [ljj with re- 
sults in general agreement with experimental results as well as more complicated theoretical 
methods (e.g., GW approach). However, at present, the application of the method is still 
limited to clusters of radius less than 1 nm. 

The study of the dielectric function via the optical transition matrix elements allows us to 
obtain information directly on absorption and photoluminescence spectra. It can also allow 
us to obtain information indirectly on the relevant radiative PL processes. The imaginary 
part of the dielectric function can be calculated by 



Ime(u;) = r- £ |e- < 0|v|5 > \ 2 5{uo - Q s ) (1) 

u s 

where e is the polarization vector of the light, v = i/h[H, r], |0 > refers to the ground state, 
\S > the excited state, and Q$ the excitation angular frequency. The excited state can be 
expanded in electron-hole pair configuration such that 



hole elec 

\s>=Y,Y, A vc\ vc > 



(2) 



The coupling coefficient A^ c can be calculated within the framework of the two-particle 
Green's function by solving the corresponding Bethe-Salpeter equation. A reasonable ap- 
proximation leads to the determination of A^ c as the positive solutions to the eigenvalue 
equation 3^ |. 



(E c - E V )A S VC + £ Ki A v , c , (n s )A s v , c , = Q S A S VC (3) 

v' ,c' 

where E c is the single-electron energy in the conduction band, E v the single-electron energy 
in the valence band, and K A C A V , C , the electron-hole interaction kernel. The most computa- 
tionally intensive part in the calculation is the evaluation of the electron-hole interaction 
kernel K AA V , C ,. This bottleneck is one of the main culprits that limits the application of the 
GW approach to only very small clusters and that of the TDLAD to clusters with the radius 
less than 1 nm. 
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Using Eq. (J2J), the optical transition matrix elements can be written as 

hole elec. 

<o|v|s>=x;£^< v M c > (4) 

v c 

The substitution of Eq. into Eq. shows that the imaginary part of the dielectric 
function depends on the transition matrix elements of electron-hole pair configuration < 
v\v\c > as well as the coupling coefficient A^ c . For clusters Si n Ge m of a given size (n + 
m=constant), the changing composition is expected to affect more on the transition matrix 
elements < v\v\c > than on the coupling coefficient A^ c because of the insensitivity of the 
EDOSs to the change in configuration (see Eq. Q and Fig. EJ). Since the main interest of 
this work is to understand the effect on the optical properties of Si n Ge m clusters of a given 
intermediate size (n + m=constant) by incorporating Ge atoms into the cluster, we therefore 
focus our attention on the effect on < v |v|c > due to the change in the composition in 
Si n Ge m clusters, particularly on how the change in < u|v|c > affects the radiative transition 
probability of the clusters as its composition varies. 

To set up a benchmark for the analysis of the trend of the the optical features of Si n Ge m 
clusters of a given size but varying compositions within the framework considered in this 
work, we have calculated the imaginary part of the dielectric function without the electron- 
hole interaction. In this situation, Eq. (jlj reduce to 

lQn 2 e 2 

Ime(cu) = — — — I < • e|c> \ 2 S(E C - E v - %uj) (5) 

* V c 

Here the sums are over all the eigenstates \v > and |c >. r is the position opera- 
tor, and V the volume of the cluster. Specifically, we calculated the average Ime(u;) 
(= E« E c |(l < v\x\c > | 2 + | < v\y\c > | 2 + | < v\z\c > \ 2 )5(E C - E v - %uj)) since the 

shapes of the relaxed Si n Ge m clusters are compact and suggesting a more-or-less isotropic 
behavior. 



Fig. presents the calculated Ime(a>) at various ratios of Si to Ge in the cases with and 
without H-passivation. It is found that the optical gaps, defined by the onset energy of the 
spectra edge in Si n Ge m clusters, are in the range of 0.4-0.5 eV and are almost independent 
of the ratio of Si to Ge. The optical spectral peaks are smooth and broad. The spectra 
display low-energy transitions and show a tail near the optical edge. The tail corresponds to 
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the defect states due to the surface distortion as shown in EDOS of Si n Ge m clusters (Fig. OJ). 
Therefore, such optical properties are not suitable for the application as optical devices. The 
optical gaps of H-terminated Si n Ge m clusters, however, are in the range of 3.3-4.1 eV and 
have a linear dependence on the ratio of Si to Ge. Unlike optical spectra of Si n Ge m clusters, 
the spectral peaks of hydrogenated Si n Ge m clusters are sharp with no tail at the edge of 
the optical spectra. Such peaks near the edge reflect the high structural symmetry (the 
tetrahedral symmetry in the present study) of the Si n Ge m H 84 clusters. The large optical 
gaps in Si n Ge m H 84 clusters indicate the possibility of optical applications. 

Since the luminescence is a result of significant overlap in electron and hole wave functions, 
and the strength of the luminescence depends on the extent of this overlap and the transition 
probability, we investigated the radiative transition probabilities of Si n Ge m and Si n Ge m H 84 
clusters. The state-to-state spontaneous transition probability (Wij) (the inverse of the 
radiative lifetime r^) of the first-order radiative process between states \i > and \j > is 
defined from the Fermi golden rule and is given by 



^ = 4 = ^ l<i|r ^" >|2 (6) 

where e and m are the electron charge and mass, respectively, c the speed of light, the 
energy difference (divided by h) between states \i > and \j >, and n the refractive index of 
Si n Ge m clusters. From ellipsometry 3^ and optical- absorption 3] experiments for porous 
Si, it seems that the refractive index decreases with increasing porosity. Since there is no 
measurement of the refractive index for Si n Ge m clusters, we choose n to be 1 in the present 
calculation. 

It is apparently from the formula of that the energy difference ufj as well as the 
dipole matrix elements < i\r\j > dominate the spontaneous transition probability. Table I 
lists the HOMO-LUMO spontaneous emission probabilities Whl (= ~~^^~\{\ < H\x\L > 
| 2 + | < H\y\L > | 2 + 1 < H\z\L > | 2 )) of Si n Ge m and Si n Ge m H 84 clusters. It can be seen that 
Whl & f° r Si„Ge m clusters are very small (~ 10~ 5 in ns^ 1 ). The corresponding radiative 
lifetime thl are quite large, i.e. 608 /is for Sin, 49 fis for Sis3Gei 8 , 71 fis for Si36Ge35, 
25 [is for Sii 8 Ge53, and 121 fis for Ge7i, respectively. But the HOMO-LUMO spontaneous 
transition probabilities Whl of hydrogenated Si n Ge m clusters are about two to three orders 
of magnitude larger than the corresponding ones of Si n Ge m clusters. Their corresponding 
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radiative lifetime thl are therefore shortened by about two to three orders, i.e., 8.1 /is for 
Si7iH 84 , 0.12 /is for Si53Gei 8 H 84 , 0.059 /is for Si36Ge35H 84 , 0.07 lis for Sii8Ge53H 84 , and 
0.098 fis for Ge 7 iH 84 , respectively. The large differences in Whl and thl between Si n Ge m 
and Si n Ge m H 84 clusters are mainly due to the energy differences of u HL between HOMO- 
LUMO states in Si n Ge m and Si n Ge m H 84 clusters. As seen in Fig. El the HOMO-LUMO 
energy differences ujhl of Si n Ge m H 84 clusters are about 10 times larger than those of Si n Ge m 
clusters. 

It is noted that the radiative lifetime t hl calculated for Si7iH 84 cluster (8.1 /xs]_Js com- 
parable to that of Si66Hg 4 cluster of a similar size (6 /zs) obtained by Hirao et al. The 
spontaneous transition probabilities Whl of Si7iH 84 cluster (0.012 xlO 4 ms" 1 ) is also con- 
sistent with the recombination rate of an excited electron-hole pair in Si crystallites (about 
10 ms 1 ) for the photon energy of 4.0 eV at 5 K calculated by Delerue et al. |6l. llOj. The 
spontaneous transition probabilities Whl is found in our calculation to be higher in pure 
Ge7iH 84 cluster (1.01 x 10~ 2 ns' 1 ) than in pure Si7iH 84 cluster (0.012 x 10~ 2 ns' 1 ). This 
result is qualitatively consistent with the results of the radiative decay rate of excitons in 
Ge quantum dots (about 0.2 x 10~ 2 ns^ 1 ) and in Si quantum dots (about 0.05 x 10~ 2 ns^ 1 ) 
with the dot radius of 10 A obtained by Takagahara et al. j^ . 

Our analysis of the radiative transition in Si n Ge m clusters indicates that the HOMO- 
LUMO spontaneous emission probability Whl (radiative life time t H l) of Si n Ge m and 
Si„Ge m H 84 clusters is very sensitive to the Ge content (m/(n + m)) in the cluster. Ta- 
ble I shows a sudden and substantial increase in Whl once Ge atoms are incorporated into 
the Si clusters. In particular for hydrogen-passivated Si n Ge m H 84 clusters, the incorporation 
of Ge atoms into the Si cluster can bring about a dramatic increase in Whl of up to two 
orders of magnitude as compared to the pure hydrogen-passivated Si cluster. An increase in 
Whl of this magnitude for Si n Ge m H 84 with n + m = 71 suggests that the incorporation of 
Ge atoms into hydrogen-passivated Si clusters may pave the way to dramatically enhance 
the optical properties of pure Si clusters. 

To shed light on the underlying physics of this dramatic change in the HOMO-LUMO 
spontaneous emission probability Whl, we carried out a detailed analysis of factors that 
might affect Whl- From Eq. ©, it can be seen that Whl is controlled by the interplay 
between the HOMO-LUMO gap ujhl and the dipole matrix-element < H\r\L >. As shown 
in Fig. |U ojhl for Si n Ge m clusters of a given size does not exhibit any sensitive dependence 
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on the Ge content. For Si n Ge m Hg4 clusters, ujhl decreases with increasing Ge content 
(m/(n + m)) and shows a linear dependence on the Ge content. However, the range of 
variation for uhl is only by a factor of ~ 1.22. On the other hand, the dependence of 
the dipole matrix element onm/(n + m) shows a more complicated pattern (see Table II). 
For example, | < H\r\L > | 2 = |(| < H\x\L > | 2 + | < H\y\L > | 2 + | < H\z\L > | 2 ) for 
Si n Ge m Hg4 clusters shows a dramatic increase of two orders of magnitude once Ge atoms 
are incorporated into the Si cluster. It peaks at the equal composition of Si and Ge atoms, 
and then reduces to a smaller but same order of magnitude for the pure Ge cluster. Since 
the range of variation for | < H\r\L > | 2 (over two orders of magnitude) far surpasses that 
for uj z hl (less than one order of magnitude), it must be the behavior of | < H\r\L > | 2 
that determines the pattern of behavior for Whl- As shown in Table I, Whl for Si n Ge m H 8 4 
clusters indeed exhibits a similar pattern as that of | < H\r\L > | 2 , namely, a drastic, almost 
two orders of magnitude, increase for small Ge content, peaking at an equal composition of 
Si and Ge atoms, reducing to a somewhat smaller value for the pure hydrogen-passivated 
Ge cluster. 

To shed light on the mechanism responsible for the dramatic increase in 
| < H\r\L > \ 2 /Whl once Ge atoms are incorporated into the hydrogen-passivated Si clus- 
ters, we express the dipole matrix element in terms of the pseudo-atomic orbitals within the 
Fireball scheme j^j, namely, 



where c^ a denotes the coefficient of expansion of the wave function A in the pseudo-atomic 
orbital a at the site i. The matrix element < ia\r\jj3 > depends only on the structural 
configuration of the cluster. Specifically 



where (p a (r — Rj) is the pseudo-atomic orbital a centered at atomic site i, Rj is the position 
vector of atomic site i, R^ = Rj — R i; and r' = r — Rj, respectively. 

Since there is no substantial structural change for hydrogen-passivated Si n Ge m H; clusters 
of a given size (n + m=constant), Eqs. (JTj) and (JHJ) then indicate that the coefficients of 



<H\r\L>=Y: < ia\r\j/3 > 



(7) 



ia,j/3 
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expansion of the HOMO and LUMO states play the most significant role in determining 
< H\r\L >. In Fig. El we plot the absolute value of the coefficient cf a for the pseudo-atomic 
orbital <p a (r — Rj) vs atomic site % for the HOMO and LUMO states of hydrogen-passivated 
Si n Ge m H 84 (n + m = 71) clusters with various Ge content. When the patterns of behavior 
for the hydrogen-passivated pure Si cluster (Si 7 iGeoH 84 ) are compared with those of the 
hydrogen-passivated pure Ge cluster (SioGe7iH 8 4), it can be seen that the HOMO state for 
the two clusters shows similar distribution. However, there is a striking difference between 
the patterns exhibited by the LUMO states of the two clusters. In the case of the pure Si 
cluster, the expansion coefficients of the p x -orbital give the major contribution to the LUMO 
states. On the other hand, two outstanding features differentiate the pattern of the LUMO 
state of the pure Ge cluster from that of the pure Si cluster: (1) The expansion coefficients 
of the s-orbital give the major contribution to the LUMO state. (2) The magnitude of the 
coefficients of the s-orbital is greater than those of the p x -orbital for the pure Si cluster 
by a factor of 2 to 3. An examination of other figures in Fig. |H1 shows that the LUMO 
state possesses the same two features as long as the cluster has a Ge content and these 
features appear mostly at the Ge sites. The overlap between a s-orbital and neighboring 
orbitals is far less restrictive as compared to that between a p a orbital with its neighboring 
orbitals. This fact, together with the more pronounced strength of the coefficients for the 
s-orbital in the LUMO state of Ge-containing clusters, leads to a stronger overlap between 
the LUMO and HOMO states of Ge-containing clusters, thus explaining the correlation 
between the appearance of these two features and the sudden and dramatic increase in the 
value of | < H\r\L > | 2 once Ge atoms are incorporated into the hydrogen-passivated Si 
clusters. This correlation also indicates that these two features are mainly responsible for 
the sudden and dramatic increase in Whl once the hydrogen-passivated Si cluster possesses 
a Ge content. 

Finally, it should be noted that the properties of the HOMO state are closely related to 
the effective mass of heavy holes while those of the LUMO state to the average effective 
mass of electrons. The mass of the heavy hole for bulk Si is reasonably close to that for 
bulk Ge. On the other hand, the average effective mass of electrons for bulk Ge is a factor 
of 2 smaller than that for bulk Si. It is therefore tempting to speculate that the behavior 
patterns of HOMO and LUMO states are related to the sameness of the mass of heavy hole 
and the difference in the effective mass of electron for Si and Ge, respectively. One may also 
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add that the fact that the effective mass of Ge is two times smaller than that of Si could 
be the contributing factor for the dominant contribution of the s-orbital in the LUMO state 
for clusters with a Ge content. 

V. CONCLUSION 

The findings of our systematic study of the structural, electronic, and optical properties 
of Si n Ge m cluster with/without hydrogen passivation include: (1) There is only a weak 
surface distortion in hydrogen-passivated clusters Si n Ge m H^ of intermediate sizes. The 
mismatch effect is then the dominating factor in determining properties of those clusters. 
(2) The HOMO-LUMO gap for Si n Ge m H; clusters of intermediate sizes opens up to 
several eVs, consistent with a blueshift in PL. It shows a linear dependence on the ratio 
n/(n + m) (n + m = constant) similar to the pattern exhibited by Sii_ x Ge x alloys [25]. (3) 
The HOMO-LUMO spontaneous emission probability Whl of Si n Ge m H; clusters is very 
sensitive to the Ge content in the cluster. Once Ge atoms are incorporated into the cluster, 
it can bring about a two-order of magnitude increase in Whl- Our analysis attributes 
this dramatic increase to the strong overlap between HOMO and LUMO states around Ge 
atoms. Our findings therefore suggest that hydrogen-passivated Si„Ge m clusters may be a 
viable candidate as components in optical devices. 
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TABLE I: HOMO-LUMO spontaneous emission probability Whl an d radiative lifetime thl of 
Si n Ge m and Si n Ge m H 8 4 clusters 



System 


W HL (ns- 1 ) 


thl (ns) 


Without hydrogenation 


Si 7 i 


0.164xl0~ 5 


6.082 xlO 5 


Si 53 Gei 8 


2.049xl0~ 5 


0.488 xlO 5 


Si36Ge35 


1.415xl(T 5 


0.706 xlO 5 


Sii 8 Ge 53 


4.078xl0~ 5 


0.245 xlO 5 


Ge 71 


0.824xl0~ 5 


1.214 xlO 5 


With hydrogenation 


Si7iH 8 4 


0.0123xl0~ 2 


81.362xl0 2 


Si53Gei 8 H 84 


0.845xl0~ 2 


1.183 xlO 2 


Si36Ge35H 8 4 


1.673xl0- 2 


0.598 xlO 2 


Sii 8 Ge53H 8 4 


1.423xl(T 2 


0.703 xlO 2 


Ge7iH 84 


1. 013X10" 2 


0.987 xlO 2 
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TABLE II: Dipole matrices | < i|r|j > | 2 = j(| < i\x\j > | 2 + | < i\y\j > | 2 + | < i\z\j > 
between HOMO-LUMO states (A 2 ) of Si n Ge m and Si„Ge m H 84 clusters. 



System | < H\r\L > 

Without hydrogenation 

Si n 0.0329 

Si 53 Gei 8 0.322 

Si36Ge 35 0.238 

Sii 8 Ge 53 0.791 

Ge n 0.181 
With hydrogenation 

Si n H 84 0.00294 

Si 53 Gei 8 H 84 0.292 

Si 36 Ge 35 H 84 0.623 

Sii 8 Ge 53 H 84 0.600 

Ge n H 84 0.464 
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FIG. 1: The stabilized structure of Si36Ge35 cluster shows a compact shape (a) and that of 
Si36Ge35H84 cluster shows a spherical-like shape with the tetrahedral symmetry in the interior 
part (b). The Si atoms are marked by the dark grey color, Ge by the black one, and H by the light 
grey one. 

FIG. 2: The average bond lengths of bsiSi, bsiGe, and boeGe in Si n Ge m clusters (a) and in 
Si n Ge m H84 clusters (c) as a function of the ratio of m/(n + m). The average angles a ^, the 
angle between two bonds (3a and (3j in Si n Ge m clusters (b) and in Si„Ge m Hg4 clusters (d) as a 
function of the ratio m/(n + m) are shown in the right side of the panel. The dotted lines in (b) 
and (d) are the tetrahedral angle of 109.47 degree in the bulk Si and Ge. 

FIG. 3: The electronic densities of states of Si n Ge m clusters (left column) and Si n Ge m H84 clusters 
(right column) at various ratio of m/(n + m). The Fermi levels are denoted by the vertical bars 
and the broaden is 0.05 eV. 

FIG. 4: The HOMO-LUMO energy gap of Si n Ge m clusters (a) and Si n Ge m Hg4 clusters (b) as a 
function of the ratio of m/(n + m). 

FIG. 5: The averaged imaginary part of the dielectric functions calculated from Eqs. (jSJ) of Si n Ge m 
clusters (left column) and Si n Ge m Hg4 clusters (right column) at various ratio of m/(n + m). The 
broaden is 0.05 eV. 
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FIG. 6: The absolute value of the coefficient cf a vs the pseudo-atomic orbital a centered at atomic 
site i for the HOMO (left column) and LUMO states (right column) of Si„Ge m H84 clusters. The 
numbers labeled at x axis denote the number of atomic site. At each atomic site, there are four 
histograms corresponding to the pseudo-atomic orbital a. The order of pseudo-atomic orbital a is 
s, p x , p y , and p z , from left to right, respectively. 
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